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Special features of surface gravity waves in deep fluid flow with constant vertical shear 



Abstract 



of velocity is studied. It is found that the mean flow velocity shear leads to non-trivial 
modification of surface gravity wave modes dispersive characteristics. Moreover, the shear 

00 ■ 

Q\ ■ induces generation of surface gravity waves by internal vortex mode perturbations. The 

o. 

, performed analytical and numerical study provides, that surface gravity waves are ef- 

"^H fectively generated by the internal perturbations at high shear rates. The generation is 



Or 

6 

> 

of view the wave generation is caused by non self-adjointness of the linear operators that 
describe the shear flow. 

PACS numbers: 92.10.Hm, 47.35. +i, 47.27.Pa 



different for the waves propagating in the different directions. Generation of surface grav- 
ity waves propagating along the main flow considerably exceeds the generation of surface 
gravity waves in the opposite direction for relatively small shear rates, whereas the later 
wave is generated more effectively for the high shear rates. From the mathematical point 
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I. INTRODUCTION 



Generation of surface gravity waves (SGW), that are the best known sea and 
oceanic waves, are naturally associated with winds. Momentum transfer from wind 
to undulating movement of the ocean which is the basic mechanism of surface waves 
generation, is investigated since Kelvin's pioneeringwork l| . Independent and inter- 
complementary theories of Fillips j^J and Miles [3, U, y, (| provide the basics of 
theoretical understanding of surface wave generation by wind. Fillips' resonant 
mechanism is responsible for excitation and initial rising of wave motion on unexcited 
surface of the fluid; Miles' mechanism — energy transfer from wind to fluid as a 
consequence of wind shear flow and surface waves interaction — is responsible for 
subsequent amplification of the waves. According to Miles' mechanism the energy 
source is wind shear flows situated outside the fluid. Other ways of SGW generation 
have been also studied such as the possibility of SGW generation by earthquakes p, 
|^ and theory of SGW generation by intra-fluid explosions In theories mentioned 
above sources of SGW generation are extrinsic for the fluid. 

The question arises as to whether intrinsic for the fluid sources (shear flows and 
vortex perturbations for example) could generate SGW. 

This question become especially interesting in view of the impressive progress 
made in the understanding of spectrally stable shear flow phenomena by hydrody- 
namic community in the past ten years. The early transient period for the pertur- 
bations has been shown to reveal rich and complicate behavior in smooth (without 
inflection point) shear flows. Particularly, it has been shown that the linear dy- 
namics of perturbations in the flows are accompanied by intense temporal energy 
exchange processes between background flow and perturbations and/or between 
different modes of perturbations. From the mathematical point of view these ef- 
fects are caused by the non-self adjointness of the linear operators in shear flows 
and are adequately described in the framework of so-called non-modal approach 



(see, e.g., jlO, |ll|, ll2(). The non-modal approach implies a change of independent 



variables from the laboratory frame to a moving frame and the study of temporal 
evolution of spatial Fourier harmonics (SFHs) of perturbations without any spectral 
expansion in time. 
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We examine the linear dynamics of surface waves and internal perturbations in 
deep fluid in the absence of wind and presence of the fluid flow with vertical shear 
of velocity. Dispersive characteristics of shear modified SGWs as well as linear 
mechanism of surface waves generation in deep fluid by internal perturbations is 
studied in detail in the framework of the non-modal approach. 

The paper is organized as follows: mathematical formalism is presented in Sec. 
II. Shear modified SGWs and their generation is analyzed in Sec. III. Applications 
of the obtained results to the concrete physical problems are discussed in Sec. IV. 
Conclusions are given in Sec. V. 



II. MATHEMATICAL FORMALISM 



Consider deep fluid with flat outer surface at z = and constant shear flow 
Uq = (Az, 0, 0) for z < 0. The shear parameter A is considered to be positive for 
the simplicity. The gravitational field is considered to be uniform g = (0,0,— g). 
Generally, four modes of perturbation: SGW, internal gravity waves, sound waves 
and vortex mode can exist in the system. To reduce the mathematical complications 
as much as possible but still keep the basic physics of our analysis, we consider fluid 
to be incompressible (which drops sound waves) and neglect the stratification effects 
(assuming that frequency of internal gravity waves is much less in comparison with 
the frequency of SGWs, i.e., considering internal gravity waves as aperiodic/vortex 
mode perturbations). In further analysis we ignore also effects of viscosity. After 
these simplifications we maintain just two modes of perturbation - SGW and vortex 
mode and write up the following differential equations for the linear dynamics of 
perturbations of velocity {u') and normalized pressure (p' = p/po): 

du' x du' 8u' z n 

— - H H = 0, (1) 

dx dy dz 

dt +AZ dx +AUz ~ dx' (2) 

d A+ Az d A = JA, (3) 

dt dx dy ' 
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with the boundary condition on the surface z = 0: 

dp' 



at 



o. 



(5) 



z=0 



Further we use the standard technique of the non-modal approach ^(|: introduc- 
tion of co-moving variables (x' — x + Azt, y' = y,z' = z, t' = t) let us to transform 
spatial inhomogeneity presented in Eqs. (jH)-([5|) to temporal one. Then, after the 
Fourier transform with respect to x' and y': 



u'(r,t) 



4vr 2 



u(k x , k y , z', t) exp [i(k x x r + k y y')] dk x dk 



the dynamic equations are reduced to: 



ik x u x + ikyUy + I — iAt k x \ u z = 0, 



(6) 



(7) 



du x 



+ Au z = -ik x p, 



(8) 



dt' 



ikyp, 



du z 
dt' 



dp 

W'~ 9Uz 



o. 



z'=0 



(9) 



(10) 



Hereafter primes of z' and t' variables will be omitted. 

From this set we readily obtain the following equation for the perturbation of 
vertical component of velocity: 



All other perturbed quantities (u x , u y and p) can be readily expressed through 
u z combining Eqs. (f7|)-([TD|): e.g., for p we have: 



d_ 



d 

— iAtk x I u z 
oz 



zAkq-zi • 



Integration of Eq. (fT2j) with respect to time yields: 

2" 



k 



d_ 

dz 



(13) 



(14) 



where F(k x , k y , z) is the constant (in time) of integration and defines the internal 
vortex mode perturbation in the flow: F(k x , k y , z) = relates to a case when the 
internal perturbation is absent. 

Fourier transformation with respect to z: 



1 

2^ 



'U'zi.kxi ky, z, t) 
F(k x , ky, z) 

reduces Eq. ()14j) to the following one: 



Uz\k x , ky, k z , t) 
F(k x , ky, k z ) 



e ikzZ dk 7 



(15) 



k (t)u z (k x , ky, k z , t) — F{k x , k y , k z ) + 4dTikC{k x , k y , t), 



where 



C 



4iirk 



d 

dz 



2iAtk x ik z J u z {k x ,ky, z,t) 



(16) 



(17) 



2=0 



Defining u z (k x , k y , k z ,t) from Eq. ()16|) . making invert Fourier transform with 
respect to k z , taking into account boundary condition \u z \ < oo at z = — oo and the 
fact that Cik z , k y ,t) does not depend on z, we obtain: 

1 f F(k x , ky, k z 



^z{kxi ky, z, t) 



2tt 



k 2 {t) 



exp(ik z z)dk z + 



+C(k x , k y , t) exp[(/c + iAtk x )z] 



;is) 



where k 2 {t) = k 2 + k 2 (t); k z {t) = k z - Atk x . 

The first term in Eq. (fTHj) relates to the vortex mode perturbation jll|, 
whereas the second term, which is exponentially decreasing with the depth, relates 
to the SFHs of shear modified surface waves. 

Substituting Eq. (fTH|) into Eq. (JT3J), and using boundary condition, Eq. (fTTj) . we 
obtain: 



d 2 C iAk r dC 



dt 2 



+ 



k dt 



+ kgC — I(k x , ky, t), 



(19) 



where 



I{k x , k y , t) 



8iA kzk 



kg 



F(k)dk z 



(20) 



x k 6 (t) k 2 (t) 

Generally, Eqs. ()19j l -(|20 p describe the dynamics of surface wave SFHs in the pres- 
ence of the internal vortical source - the term I(k x , k y , t) is the result of an interplay 
of the mean flow shear and the internal vortical perturbations and couples the later 
perturbation with the surface one. So, there is no coupling between these pertur- 
bations in the absence of the shear. Indeed, if there are no surface perturbations 
initially [u z {k x , k y , z = 0, t = 0) = 0], then from Eqs. (fTfij) and (|20j). one readily ob- 
tains that I(k x , k y , t) ~ u z {k x , k y , z = 0, t = 0) at A = 0, i.e., I(k x , k y , t) = 0. Thus, 
if there is no the source in shearless flow initially, it does not appear afterward. 



III. SGWS AND THEIR GENERATION IN SHEAR FLOW 

One can see from Eqs. (|T9|) and (J20|) . that there are two main effects of the shear: 
firstly, the second term on the left hand side of Eq. indicates that the velocity 
shear affects the frequencies of SGWs. Secondly, the source term I(k x , k y ,t) caused 
by the internal perturbations, couples the internal and surface perturbations and 
results the emergence/generation of SGW in the flow. Our further attempts are 
focused on the study of these effects. 



A. Shear modified SGWs 



In this subsection we study shear induced modifications of properties of SGWs. 
For this purposes we assume that initially there were no vortex mode perturbations: 
F(k x , k y , k x ) = 0. Consequently, I(k x , k y , t) = [see Eqs. (|'2L)|)]. and Eq. (jTUj) reduces 
to the homogeneous one, with the solution: 

C h {k x , k y ,t) = Ci{k x , k y ) exp(-ifM) + C 2 {k x , k y ) exp(-iQ 2 t), (21) 

where Ci^{k x , k y ) are determined by initial conditions and 

n - = ± ^W - f = ^ ( ± /^ - 4) (22) 
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represents shear modified frequencies of SFH of SGW propagating in opposite direc- 
tions and S = AI(Akq) 1 / 2 is the dimensionless shear rate. This equation provides, 

rnn 

that in contrast with acoustic and magnetohydrodynamic wave modes |14L 1151 116|. 
the presence of the shear does not lead to the time variability of the frequency. 
However, velocity shear leads to the non-trivial modification of the frequencies and 

17, Q|. Indeed, for the value of the phase 



consequently phase velocities of SFH 
velocity, Eq. (J22J) provides: 



cos 2 — Scos<f)) , (23) 
V k ^ ' 

where = arccos(A; x //c). 

The phase velocity is isotropic in the shearless limit (S = 0), while depends on 
in the shear flow. The anisotropy increases with the shear rate. The value of the 



phase velocity is minimal at = 0: V^ m = y g/k(yl + S 2 — S) and is maximal 
at <fi = n: V™ x = y g/k(y/l + S 2 + S). Suppose SGW is emitted by point source 



situated on the surface at x — y — 0. From Eq. (|23j) it follows that the propagation 
of the leading wave crest is described by: 

r(S } M = V ph (S,(f>)t = 

= Jl [yi + S 2 cos 2 - S cos 0) t. (24) 

Figure Q shows the leading wave crest of the SGW for three different moments of 
time ti, t2, ts, when t2 = 2t\, t$ = 3ti, that are circular but not concentric. 



B. Generation of SGWs by internal vortices 

Let us initially analyze the source term I(k x ,k y ,t) that is determined by 
Fiykxi kyj k^) • Assume F(^kx^ ky^ k% ) is a localized function in the wavenumber space, 
with a center of localization at k Q = (k x o, k y o, k z o). Note that the first multiplier in 
the integrand of Eq. (|20|) reaches its maximum when k z — Ak x t = 0. Consequently, 
the maximum of the integral is in the vicinity of time t — t* = k z o/(Ak x o). Equa- 
tion ()20)1 provides, that generally I(k x , k y , t) tends to zero in both limits t — > ±00. 
Actually, it exists some time interval 2 At around i* where the source term differs 
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from zero. The value of At depends on the degree of localization of internal per- 
turbation, i.e., of F(k x , k y , k z ) in wavenumber space. (The source localization is 
demonstrated below on a specific example.) Thus, in the case of a localize source, 
the coupling between surface (gravity wave) and internal (vortex mode) perturba- 
tions takes place in some time interval 2At around and at \t — t*\ > At this 
perturbations can be considered separately. 

The general solution of the inhomogeneous equation, Eq. (j!9|) . is the sum of the 
general solution of corresponding homogenous equation and a partial solution of the 
equation: 

C(k x , ky, t) = Ch{k x , ky, i) + Ci{k x , k y , t). (25) 

The general solution Ch{k x , k y , t) is presented by Eq. (|2*T]) . whereas a partial solution 
of Eq. (JH is: 

1 ft 

Ci = — — exp (— ittit) / I(k x , k y , t')exp (iVtit')dt' 

i r* 

- — — exp {-iVL 2 t) I I(k x , k y , t') exp {iVt 2 t')dt' , (26) 
2"o J to 



where 



kg + 



A 2 k 2 x 
4k 2 



k 2 

'hjJl + S 2 ^. 



(27) 



Assume that the coupling between the surface and internal modes can be ne- 
glected at the initial moment of time to, i.e., t < t* — At. After passing through 
the coupling time interval, for any t > t f = t* + At the modes become independent 
again. However, during the time interval [^o^/L internal vortices generate SGWs 
with frequencies Q± and fl 2 [see Eq. f!22j) ]. As it follows from Eqs. (f2~Tj) . (|2*K|) . and 
(}2l)j) . if initially there are no SGWs (Ci j2 = 0), the generated SFH amplitudes (Qi^) 
are: 



Qi(k x , ky 



2ttr 



I(k x , k y , t') exp (iflit')dt' 



to 



(28) 



Q 2 {k x , ky 



2Qr 



I(k x , k y , t') exp (iQ 2 t')dt' 



to 



(29) 



One can replace the limits of integration by ±oo. After integration in time this 



S 



yields: 



Q 



1,2 



7tkx 
k 3 



y4. Jvq 
2Q k 



F[k x , k y , k z ) exp 



2ir 2 k x 
k 3 



y4. kx 



(ilo =F A/2)k 

^n oT A/2)k z 
(n oT A/2)k 



dk 7 



x F ( k x , ky, 



T A/2 



(30) 



Note that the last multipliers in Eq. ()3()|) are proportional to the vorticity of 
initial perturbations at Z\ t % = — (f2 =F A/2)/(Ak x ) respectively. The second mul- 



tipliers indicates, that at small shear rates (S = A/y4kg 1), the amplitudes 
of generated SGWs are exponentially small with respect to the large parameter 
1/5". Equation (}3T)j) also indicates that for a fixed k x , the generation has maximal 
effectiveness in two dimensional case (k y = 0). 

Let us describe the dynamic picture for a specific example, when pure internal 
vortex mode perturbation (without any admix of surface waves) is imposed in the 
flow initially. For simplicity we consider two dimensional problem, when d/dy = 0. 
The vertical velocity of the imposed perturbation is given by: 



u z (x,z,t ) = z i r](-z)exp 



[(z + z ) cos + x sin . 



x exp 



[(z + zo) sin — x cos ( 



(31) 



where n(z) is Heaviside function, (0, —z Q ) is the center of the localization, L 12 char- 
acterize vertical and horizontal scales respectively and <fi is the slope of the pertur- 
bation. 

Numerical solution of the problem was performed as follows: Fourier transform 
of Eq. (|3*T|) with respect to x variable allows to determine F(k x , z) through Eq. ()14|) . 
Another Fourier transform with respect to z yields F(k x ,k z ). Then, the source 
function I(k x ,t) is found by Eq. (|2(Jj) . Thus, the solution of the problem for a fixed 
k x reduces to the numerical solution of the inhomogeneous equation, Eq. (fTU|) . with 
known I(k x , t). 

Dependence of the source function I(k x ,t) on t at L\ = 1, Li = 7, (f> — 7r/18, 
k x — 1, z — 2 for two different values of the shear rate 5* = 0.08 (dashed line) and 



S = 0.32 (solid line) is presented in Fig. El As it was mentioned above, the source 
term is localized function and considerably differs from zero only in the interval 
t E (20, 40) for S = 0.08 and t E (5, 10) for S = 0.32. 

For analysis of the wave generation effectiveness it is useful to introduce the 
generation coefficients, that characterize the ratio of generated wave energy density 
and the maximum energy density of the initial vortex mode perturbations for a fixed 
value of k x . Taking into account that the maximum energy density of the vortex 
mode perturbations is: 

1 f 00 

Ev = Wj \ F ^ z )\ 2dz ^ ( 32 ) 
and the energy density of the generated waves: 

E w ip, = -r-Qi^i^x), (33) 
we define the non-dimensional generation coefficients as: 

Gi,2 = ( joo^ \F(C, z)\ 2 dz J ■ W 

Figure H3 represents generation coefficients G\ (dashed line) and G2 (solid line) 
vs shear rate S at Li = 1, L2 = 7, <f> = tt/18, k x = 1 and zq = 2. As it can 
be seen, at small values of the shear rate, generation of SGW with frequency Q± 
(i.e., propagating along x axis) considerably exceeds the generation of SGW with 
frequency ^2 (he., propagating against x axis), whereas the later wave is generated 
more effectively at S > 0.15. 

The wave generation is well traced in Figs. 0] and where the temporal evolu- 
tion of the vertical component of velocity perturbation at the surface obtained by 
numerical solution of Eqs. (fTT)|) - (|2T)|) is presented for S = 0.32 and S = 0.08 respec- 
tively. The other parameters are the same as in Fig. |21 Purely internal vortex mode 
perturbation is imposed in the equations initially. The generation takes place in 
the time interval where I(k z ,t) noticeably differs from zero. Afterward, just (two) 
waves with different frequencies and amplitudes exist. At S = 0.32, presented in 
Fig. HI the generation takes place in the time interval t E (5,10). Besides, SGW 
propagating against x axis is mainly generated. In contrast to this, at S = 0.08, 
presented in Fig. the generation of SGW propagating along x axis is dominated. 
These numerical results are in agreement to analytical ones [see Eq. f!30|) and Fig. |3] . 
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IV. DISCUSSION 



In the previous sections simplified model was considered that allowed us to sim- 
plify mathematical description and study shear induced effects in the 'pure' form. 
For instance, the influence of the viscosity was ignored and the density ratio p a j p 
of the fluids above and below of the surface z = was assumed to be zero. The 
later assumption allows to ignore all the dynamical processes in the upper fluid. 
On the other hand, it is well known that in the case of ocean waves the wind is the 
most important and powerful source of the waves. In this section we discuss possible 
applications of the studied linear effects to the concrete physical applications. 



A. Ocean waves 

It is well known \A that the wind is main source of ocean SGWs. 

In the context of future discussion the papers of Chalikov group [23j, |24| should 
also be noted, where the influence of small scale turbulence in the air on wave 
growth was studied in detail. At present there exists well developed theory both 
SGW generation and nonlinear evolution that is mainly confirmed by experiments 
as well as numerical simulations (for recent review see, e.g., After development 

of wind driven instability nonlinear 4-wave resonant interactions transfer the wave 
energy to smaller scales. Existing theory predicts that for relatively small frequencies 
Zakharov-Philonenko j3] spectrum E(u) ~ uj~ 4 of SGW fluctuations (sometimes 
called Toba's spectrum) should be observed (in this context see also j^ll]). whereas for 
relatively high wave numbers nonlinearity becomes strong and Phillips's spectrum 
E(cj) ~ co>~ 5 of the wave turbulence should develop. Existing observations confirm 
this predictions and provide that in the range u p /3 < uj < 3u p , where uo p is the 
peak frequency, Zakharov-Philonenko spectrum is usuahy observed. For „ > 3, p 
the spectrum becomes very close to Phillips's one [19||- The properties of the wave 
spectrum in the very short wavelength region, as well as dynamics of SGW turbulent 
fluctuations dissipation is much more unclear 2^ |. 

In the case of ocean waves presented linear mechanism of SGW generation can 
have important contribution to the balance of small scale SGW fluctuations. Indeed, 
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characteristic length scale of the turbulence at the ocean surface is much smaller than 
in the air. Namely, the characteristic length and velocity scales are w* ~ 1 cm/ sec 
and I ~ 1 cm presence of the wind strong 

velocity shear A ~ 10 sec~ l is presented in so called 'buffer layer' j3] of the water, 
with the thickness l\ ~ (20 — 100)Zo, where Iq ~ v ju* is the dissipation length scale 
and v is kinematic viscosity of the water. Simple estimates yields l\ ~ (0.5—1) cm. 
Presented linear mechanism implies that vortical perturbations generate SGWs with 
the same length scale. Therefore in the case of ocean waves internal vortex mode 
perturbations should effectively generate small scale SGWs - with the wavelength 
just above the capillary length scale A c = 0.39 cm j^. In this context the study 
of the influence of capillary effects on the processes discussed above seems to be 
interesting. Analysis of this problem will be presented elsewhere. 



B. Interfacial gravity waves 

In the analysis presented in Sees. II and III density ratio p a /Po of the fluids 
above and below of the surface z = was assumed to be zero. Obtained results can 
be readily generalized in the case of interfacial GWs. If the density of upper and 
lower fluids are p\ and p 2 and the shear rates are A\ and A2 respectively, then shear 
modified dispersion of interfacial GWs is given by the same expression (122)1 with g 
and A replaced by g* and A*, where 

P2-P1 , Ap2 - Axpi 

9* = 9 ; , A* = ; • (35) 

Pi + Pi P2 + Pi 

From this equation it follows, that the influence of shear both on the wave dispersion 
and coupling with internal vortex perturbations, that is determined by dimensionless 
parameter 

c - A * q 1 ~ P1M/P2A2 . , 

v 1 - pVpI 

is much more notable when the fluids have comparable densities if p\A x is not very 
close to p2^2- Therefore, described shear induced effects usually should have much 
more influence on the dynamics of interfacial waves, then on the ocean waves. 
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V. SUMMARY 



Let us summarize the main features of the linear dynamics of surface gravity 
weaves in a simplified deep fluid (at z < 0) flow with vertical shear of the mean 
velocity Uq = (Az, 0,0). The simplification lies in the neglecting of the fluid com- 
pressibility and stratification, in other words, in the consideration of the system 
containing just two modes of perturbation: surface gravity wave mode and internal 
vortex mode. Special features of SGW in the system are the following: 

The mean flow velocity shear causes the non-trivial modification of the frequencies 
and phase velocities of SGWs. The frequencies are defined by Eq. ()22p. The phase 
velocity becomes anisotropic [see Eq. (J23j) and Fig. Q: its value is minimal for 
SFH propagating along x axis [V™ in = y / ^/fc(v / TT52 - S)} and maximal for SFH 
propagating against x axis [V^ m = \j gjki^JX + S 2 + S)]. 

The mean flow velocity shear leads to the appearance of the intrinsic (to the 
fluid) source of SGW generation via coupling the wave with the internal vortex 
mode perturbations - the coupling results the emergency /generation of SGWs by 
internal vortex mode perturbations at S > 0.05. The generation is different for the 
waves propagating in the different directions [see Eq. OOJ - Generation of SGW with 
frequency Qi considerably exceeds the generation of SGW with frequency Q2 for 
relatively small shear rates S, whereas the later wave is generated more effectively 
for the high shear rates (S > 0.15). 
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FIG. 1: Shear induced anisotropy of SGW propagation: The leading wave crest at 
three different moments of time t±, ti, t%, when ti = 2t±, t$ = 3ii, that are circular, 
but not concentric. A point source of the SGW is located at x = y = 0. 
FIG. 2: I(k x ,ky,t) vs time at S = 0.32 (solid line) and S = 0.08 (dashed line) 
k x = 1, L\ = 1, L2 = 7, zq = 2 and <fi = -k /IS. 

FIG. 3: Generation coefficients G\ (dashed line) and G2 (solid line) vs. shear rate S 
at k x = 1, L\ = 1, L2 = 7, zq = 2 and <f> = tt/18. 

FIG. 4: u z (k x ,t) vs. time at 5 = 0.32, k x = 1, L ± = 1, L 2 = 7, z = 2 and 4> = vr/18. 
FIG. 5: u z (k x ,t) vs. time at S 1 = 0.08, k x = 1, L± = 1, L2 = 7, zq = 2 and = 7r/18. 
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